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Abstract 

We provide evidence of a superfluid-insulator transition (SIT) of magnons in a quasi-onc-dimensional 
quantum ferrimagnet with isotropic competing antiferromagnetic spin interactions. This SIT occurs between 
two distinct ferrimagnctic phases due to the frustration-induced closing of the gap to a magnon excitation. It 
thus causes a coherent superposition of singlet and triplet states at lattice unit cells and a power-law decay on 
the staggered spin correlation function along the transverse direction to the spontaneous magnetization. A 
hard-core boson map suggests that asymptotically close to the SIT the magnons attain the Tonks-Girardeau 
limit. The quantized nature of the condensed singlets is observed before a first-order transition to a singlet 
magnetic spiral phase accompanied by critical antiferromagnetic ordering. In the limit of strong frustration, 
the system undergoes a decoupling transition to an isolated gapped two-leg ladder and a critical single linear 
chain. 
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I. INTRODUCTION 



Recently, several experimental and theoretical studies indicate that, under very special condi- 



tions, magnons 



Q.y. 



and polaritons ^\ undergo Bose-Einstein condensation (BEC) in two- and 



three-dimensional materials. In magnetic systems, BEC of magnons can be driven by an applied 
magnetic field (h) (Ref. [ll), by varying the external pressure [2I, or by microwave purnping [s] In 
ID gapped antiferromagnets, e. g., spin-1 chains Q] and single spin-1/2 two-leg ladders [g], the gap 
to the magnon excitation closes at a critical value {he) of the field and the magnetization increases 
as {h — hc)^/'^ . Although, stricto sensu, there is no BEC of magnons in these ID systems, it is very 
appealing to describe the transition in terms of the condensation of the uniform component of the 
magnetization along the applied field Q]. In fact, rigorous results tI on low dimensional {D < 2) 
uniform interacting boson systems preclude the occurrence of BEC in finite temperature (T). In 
2D systems phase fluctuations have mainly a thermal origin, so that only the T = condensate 
survives, with superfluid behavior persisting up to the Kosterlitz-Thouless temperature. In con- 
trast, in ID boson systems phase fluctuations have a quantum origin and there is no BEC, even at 
T = 0, but superfluidity is expected Q]- However, in finite systems the scenario is more complex, 
since in real confined systems [3, one may be dealing with metastable states. 

In this work we introduce an isotropic Heisenberg spin Hamiltonian with two competing anti- 
ferromagnetic (AF) exchange couplings [Ji(= 1) and J] exhibiting a continuous quantum phase 
transition at a critical value Jd which, we argue, is a superfluid-insulator transition (SIT) of 
magnons associated with the creation of a coherent superposition of singlet and triplet states at 
lattice unit cells. For J = 0, the model shares its phenomenology and unit cell topology with 



quasi-one-dimensiona! 



ferrimagnetic compounds Q], such as the line^ trimer clusters present in 



copper phosphates [10(], and the organic ferrimagnet PNNBNO (Ref. On the theoretical side, 

several features of the ferrimagnetic phase have been studied through Hubbard [l^, t — J (Ref. 
[1^ ) and Heisenberg [l^ models, including magnetic excitations [iS, [3| and the occurrence of 



new phases induced by hole doping of the electronic band 



171 ]. Also, the physical properties of the 



compound CU3 (003)2 (0H)2 were successfully explained [18|] by the distorted diamond chain model 
[3], which is a system with three spin 1/2 magnetic sites per unit cell and coupling parameters 
such that the ferrimagnetic state is frustrated. 

Numerical results have been obtained for finite clusters through Density Matrix Renormalization 
Group (DMRG) (Refs. 20|, |2l[) using open boundary conditions and exact diagonalization (ED) 



using periodic boundary conditions and Lanczos algorithm. 
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FIG. 1: (a) Illustration of the A and B sublattices (circles) and AF spin couplings which favor (full lines) 
and destabilize (dashed lines) the LM ferrimagnctic GS: Ji(= 1) and J, respectively, (b) Illustration of the 
LM ferrimagnetic GS. (c) Results (see text) for Sg/SLM] dashed and dotted hues are guides to the eye. 



The paper is organized as follows: in Sec. [TI]we introduce the model Hamiltonian and analyze 
the magnetic correlations of the competing phases close to J = Jd- In Sec. IIIII we define a 
hard-core boson model (HCB model), which is used to describe the main characteristics of the 
magnon SIT at J = Jd, in particular, the Tonks- Girardeau hmit. Further, in Sec. IIVI we discuss 
the singlet magnetic spiral phase accompanied by critical antiferromagnetic ordering, which sets 
in after a first-order transition at J = J^, as well as the decoupling transition, at J = Jc2, to an 
isolated gapped two-leg ladder and a critical single linear chain. Finally, a summary of the results 
is presented in Sec. |Vl 

II. MODEL HAMILTONIAN AND ORDERED PHASES 

The model Hamiltonian reads: 

1=1 a=l,2 I 
+Bii • B2Z + ^ Ba,l ■ Bq,_;+i), (1) 

Q = l,2 
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FIG. 2: DMRG results for the magnetic structure factor, Fx{q), at q = and q = tt for {&) X = A and (b) 
X = B (Bi or B2) spins in a chain with Nc — 33; dashed hues are guides to the eye. (c) lUustration of the 
F2 phase. 



as sketched in Fig. ^a). In Eq. ([T]), A;, Bi; and denote spin 1/2 operators at sites A;, Bi; 
and B2/ of the unit cell /, respectively, and Nc is the number of unit cells. For J = the model 
(named AB2 chain or diagonal ladder) is bipartite and the Lieb-Mattis (LM) theorem 23] predicts 
a ground state (GS) total spin 



\Na-Nb\ 



Slm, 



(2) 



2 2 

where Na (Nb) is the number of A (Bi and B2) sites. The GS spin pattern is represented in Fig. 
[U^b). In Fig. [T^c) we report data for Sg/SLM as a function of J using DMRG {Nc = 33) and ED 
(Nc = 10). Although the LM theorem is not applicable for J 7^ 0, the ferrimagnetic phase {Fl 
phase) is robust up to J w 0.342 = Jd (Ref. [23|]), beyond which Sg steadily decreases {F2 phase) 
before a first order transition to a phase with Sg = (apart from finite size effects) at J ~ 0.445. 
In order to characterize the F2 phase, we have calculated the magnetic structure factor, 

Fx{q)=^Cx{l)e"^\ (3) 
/ 

with q = 2TTn/{Nc — 1), n = 0,1,. ..,Nc — 1, where Cx{l) is the two-point correlation function 
between spins separated by / unit cells at sites X = A,Bi and B2. We first noticed that the A 
spins remain ferromagnetically ordered as the critical point Jd = 0.342 is crossed, although the 
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FIG. 3: DMRG results for the square of the longitudinal (L) and transverse (T) order parameters at the 
spin se ctor 5^ = Sg and J = 0.395 for (a) A and (b) B (Bi or B2) spins (full lines are polynomial fittings), 
(c) DMRG results for the transverse staggered correlation function Cbttt{1)- 

magnitude of the peak at g = decreases for J > Jd, as displayed in Fig. [2]^a), while no peak 
is observed at g = vr. The Bj (i = 1 or 2) spins also remain ferromagnetically ordered (peak at 
q = 0), with similar J-dependence, as shown in Fig. ^h). However, an extra peak at q = tt 
develops after the transition, which indicates the occurrence of a period-2 modulation in the spin 
pattern for J > Jd. Further, the average value of the correlation function (Bi^ • B2/), which 
amounts to ~ 0.25 (triplet state) in the Fl phase, steadily decreases after the transition at Jd- 
These findings suggest that the F2 phase would display a canted configuration, as illustrated in Fig. 
[2]^c). However, to check whether these features are robust in the thermodynamic limit, we have 
studied the finite size scaling behavior of the transverse (T) and longitudinal (L) order parameters 
in the F2 phase: 



for q = (uniform component) and q = tt (staggered component), in the subspace of maximum to- 
tal spin z-component {S^ = Sg). The correlations are studied at J = 0.395, for which Sg = Slm/2, 
and the results are shown in Fig. [3l We confirmed that in the (extrapolated) thermodynamic 



m 




Fx(L.T) {q) 



(4) 
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limit the spins at sites A and B are ferromagnetically ordered, as indicated by m\j^{q = 0) 7^ in 
Figs. EJa) and (b). Further, since the A and B net magnetizations are oppositely oriented, the F2 
phase is ferrimagnetic. The values of 'm\j^{q = vr), m'^j^{q = vr) and m?^rp{q = 0) nullifies linearly 
with system size, which evidences short-range correlations. On the other hand, the best fitting to 
the data for m^j.(g = vr) presents a nonlinear dependence with the inverse of the system size and 
also nullifies in the thermodynamic limit. This behavior indicates that the staggered correlation 
function of the spins at sites B along the transverse direction to the spontaneous magnetization, 
Cbt-k{1)-, exhibits a power-law decay, as explicitly confirmed in Fig. [3|^c). We thus conclude that 
for Jci < J < 0.445 the GS is also ferrimagnetic but with critical correlations along the transverse 
direction to the spontaneous magnetization (F2 phase). 

Next we focus on the effect of J on the magnetic excitations. For J = the Hamiltonian 
exhibits three magnon modes {15I . 1^. One is AF, i. e., the spin is raised by one unit with respect 
to the GS total spin, while the other two are ferromagnetic, associated with the lowering of the 
GS total spin by one unit. The AF gapped dispersive mode is responsible for a quantized plateau 
in the magnetization curve as function of h and should also exhibit condensation, as suggested 
by the numerical data in Ref. [3]. One of the ferromagnetic magnons is the gapless dispersive 
Goldstone mode, while the other is a flat mode and is the relevant excitation for the transition 
at J = Jci- To understand some nontrivial features of this excitation, we must comment on the 
symmetry properties of the model. For J = the Hamiltonian is invariant under the exchange 
of the B sites at the same cell. This symmetry implies that spins at these B sites can be found 
only in singlet or triplet states (mutually exclusive possibilities); in the GS only triplets are found. 
The relevant magnon is a localized gapped mode which induces the formation of a singlet pair in 
one cell, as illustrated in Fig. [D^a). For J 7^ 0, this local symmetry is explicitly broken and the 
spins at these B sites can be found in a coherent superposition of singlet and triplet states. In Fig. 
lUb), data using ED for the magnon band (g = 2Tin/Nc, ?i = 0, 1, A'^c — 1) is displayed for various 
values of J before the transition point. For J = the band is flat with a gap Aq ~ 1.0004. By 
increasing J, the bandwidth increases and the gap to the GS lowers, closing at the wave vector 
(7 = vr at the transition point. 
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FIG. 4: (a) Illustration of the relevant magnon excitation for J = 0: ellipse indicates a local singlet state, 
(b) Magnon band for J = 0.00,0.05,0.10,0.15,0.20,0.25,0.30 and 0.34, from top to bottom, (c) Ns as a 
function of J. Full lines are the HCB model predictions in the TG limit and dashed lines are guides to the 
eye. 
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III. HARD-CORE BOSON MODEL, SUPERFLUID-INSULATOR TRANSITION AND 
THE TONKS-GIRARDEAU LIMIT 



The GS total number of singlets is given by 



Ns = J2^m), (5) 

with singlet density (rji) = (sjsi), where 



1=1 



is the creation operator of a singlet pair at cell / and Bj^ ^ is the creation operator of an electron 
with spin a at the Bj (i = 1, 2) site of cell I. In fact, it is easy to show that 

{m) = \-{Bn-B2i), (7) 

so Ns = for J = 0. In Fig. \MyC) we observe that Ns starts to increase in steps of unity after 
J = Jci, indicating the quantized nature of the condensing singlets. 

We now examine the nature of the quantum critical point at J = Jd- For this purpose we split 
the Hamiltonian of Eq. ([1]) in three terms: the first favors ferrimagnetism, 

Hab = Y.Ai-{Si + Si.i), (8) 

where Si = Bi; + B21; the second one favors AF ordering between A spins, i. e., 

Ha = JY,Ai-Ai+^, (9) 

I 

and shall play no significant role in our analysis; the last term, also unfavorable to ferrimagnetism, 
is a two-leg ladder Hamiltonian connecting spins at sites Bi and B2 (discarding a constant factor) 



24 



= ^ ( E + E • S'+i + E • Di+i ) , (10) 
\ I I I / 



where D; = Bi; — ^21- We represent the Hamiltonian in a basis with two states for each pair 
Bii and B21: the singlet and the triplet component in the magnetization direction. In addition, 
we define the vacuum of the HCB model as the state with this triplet component in each cell. 
We now study the GS energy when a number Ns of singlet pairs is added to the vacuum. For 
J = 0, the energy cost of a singlet pair is Aq (the gap to the flat mode); thus, for Ns singlets, the 



contribution from H^b is NsAq. The first term in Hb is diagonal and will add a factor of —JNs; 
the second causes a repulsion between singlets and adds also an extra factor of —JNs; finally, 
the last term in Hb introduces the singlet itinerancy. Grouping these contributions, we arise to a 
model of hard-core bosons with nearest-neighbor repulsion: 



Hs = (Ao - 2J)Ns + ^Y1 'm+i + ^ Y.^s\si+i + h.c). 



(11) 



/ I 

We remark that the hard-core boson interaction is implied by the algebra of the singlet operators: 



1; and 

for / 7^ m. 



Before the transition, the single magnon dispersion relation. 



oJq{J) = Ao — 2 J + J cos 



(12) 
(13) 



(14) 



agrees well with the numerical data for g ~ vr, as can be seen in Fig. [D^b). The resulting critical 
point: ujq=-n{Jci,s) = 0, i. e., 



Jci,s = 0.333, 



(15) 



is in excellent agreement with the numerical prediction Jd = 0.342. Moreover, the closing of the 
magnon gap is also in excellent agreement with the prediction 



Aj = u:{q = vr) = 3(Jci,s - J) 



(16) 



and with the expected linear vanishing of the Mott gap 26|]: zv = \, where z = 2 and u = \/2 are 
the correlation length and dynamic critical exponents, respectively (see below). 

After the transition and in the highly diluted limit (-^ = r/ — > 0), the energy of Ns hard-core 
bosons in ID is well approximated by the energy of Ns free spinless fermions [5]. Through this 
map, the energy density reads: 

Egs{J) ^ f"' 

-kp 



£gs{J) 



3(Jci,s - J)il + 



dk, , ^ 



6 



(17) 
(18) 



where kp = ttt] and 



ek{J) - fJ-F = ^k+-n{J) ~ -3(J - Jci,s) + 



(19) 
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FIG. 5: (a) Ground state energy, EqSi relative to the energy of the LM state, (b) Luttingcr Liquid exponent, 
as function of r}. Full lines are the HCB model predictions in the TG limit and dashed lines are guides 
to the eye. 



2£| 



Notice that the Fermi chemical potential satisfies the Tonks- Girardeau (TG) limit [3, [27I . 
(ID Bose gas of impenetrable particles), corresponding to an infinitely high repulsive potential in 
the Lieb-Liniger solution [3] of the 5-function ID Bose gas: 

7r2 Jr?2 



= ep{J) 



(20) 



where J ^ is the fermion mass, h = 1, and i] is the density of singlets for J > Jci,s derived from 
the equilibrium condition dr/ScsiJ) = 0: 

_ V6(^ - Jci,s) 



-, 7?^0, 



(21) 



much in analogy with the ID field-induced transition. Further, in Fig. Hlld) we display the good 
agreement between the numerical estimate for the density (rji) of a two (four) particle state, Ns = 2 
{Ns = 4), in an open system and the HCB model in a continuum space given by 24 1 



N, 



-^sin^(A:„/), 



(22) 



n=l 



with kn = 1, -^^TTf- Also, as shown in Figs. [U^c), Hl^c) and[5]^a), the HCB model predictions for 



So 



Sl 



M 



l-2ry, 



(23) 



r] and Eos (J), respectively, are very close to the numerical data for J > Jci,s- 



On the other hand, using the Luttinger liquid description 



3C, 



3ll | for our highly diluted HCB 



10 



model we have the following general relations for the sound velocity c and the compressibility k: 

-kJt] 



1 



K 

ire 



(24) 
(25) 



where K is the Luttinger parameter governing the decay of the correlation functions. However, 
since 



1 (j'SoS 2 J 



(26) 

it implies that K = 1; thus c = vrJ?] = Jkp, in accord with the TG limit [3, Further, taking 
r] as the order parameter of the SIT, Eq. (j2ip implies (3 = 1/2, while r/^K diverges with a critical 
exponent = 7 = 1/2, in agreement with the scaling and hyperscaling relations 261 ): 



a + 2/3 + 7 = 2, 

2 — a = i^((i + 

respectively, assuring that the SIT is in the free spinless gas universality class 
In an interacting Bose gas 



(27) 
(28) 



32|. 



331], K = 1/2 is the separatrix between systems dominated by 



uctuations, K < 1/2, 



34| have succeeded in 



superfluid fluctuations, K > 1/2, from those dominated by charge density fli 
(in our magnetic model spin fluctuations prevail). Affleck and collaborators 
taking into account corrections from interactions between pairs of dilute magnons parametrized by 
a scattering length, a, thus implying that 

K = l-2am + 0{m^), 



where m is the field-induced magnetization for the S = 1 chain, with a 
increase of K with m was confirmed by numerical calculations 



(29) 

-2. The predicted 



34i |. This parametrization can also 



be implemented in our problem. In fact, in Fig. [5^b) we show that K = 1 — At], with a ~ 2, 
fits quite well the data for the Luttinger liquid parameter in the highly diluted regime. K was 
calculated using DMRG and assuming 



^ 00 

I2K 



(30) 



IV. SPIRAL CORRELATIONS, WEAKLY COUPLED AF CHAINS, AND LADDER- 
CHAIN DECOUPLING 



We now turn our attention to the transition point Jt ~ 0.445, which marks the onset of a singlet 
phase, as can be seen in Fig. 1(c), characterized by non-quantized values of Ns, as shown in Fig. 

11 



4(c). On the other hand, from the Hamiltonians in Eqs. (8)-(10), we can infer that for J >> 1 the 
system should decompose into a linear chain (A sites) and an isotropic two-leg ladder system (Bi 
and B2 sites); see Fig. 1(b). The linear chain is known to be gapless with critical spin correlations 
(power-law decay), while the two-leg ladder is gapped with exponentially decaying correlations. In 
what follows we discuss the complex phase diagram in the region J > Jf. 

Initially, we display in Fig. 6 the magnetic structure factors F^^q) and Fsiiq), with i = 1 or 2, 
as well as Fs{q), which is associated to the magnetic structure of the composite spin Si = Bi; -|-B2z. 
In Fig. 6(a) we see that FA^q) peaks at q = for J = 0.44, i. e., the system remains in the F2 
phase and the A spins are ferromagnetically ordered. For J = 0.45, a sharp peak in a spiral wave- 
vector qmax is observed. The peak broadens and qmax increases with increasing J. For J = 0.56 we 
notice the emergence of a commensurate AF peak, coexisting with the spiral one, particularly for 
J > 0.60, as seen in Figs. 6(a) and 6(b). On the other hand, we observe in Fig. 6(c) the presence of 
two peaks in Fs^iq) for J = 0.44: the q = peak associated with the ferromagnetic ordering of the 
Bi sites in the F2 phase, and the q = tt peak related to the critical staggered transverse correlation 
at the same phase. Likewise, for J = 0.45, a spiral peak is observed at the same wave-vector qmax 
of FA{q)- Further, notice in Fig. [6)[d) that the magnitude of the AF peak drops in the interval 
0.96 < J < 1.00. 

In order to develop a physical meaning of the above referred data, we first point out that the 
coupling between spins at A and B sites occurs through the composition Si = Bi; -|- B2/, as can 
be seen in Eq. (8). Further, as the singlet component of Si is magnetically inert, only its triplet 
components affect the magnetic ordering at the A sites. In fact, as shown in Figs. 6(e) and 6(f), 
short-range spiral ordering is observed in the magnetic structure of Si up to J ~ 1.00. However, 
since the peak is weak and broad for J > 0.6, its feature is overcomed by the AF one in the data 
of Figs. 6 [(a)-(d)]. In the sequence, we focus on the AF ordering observed for J > 0.6 and study 
how the system approaches the ladder-chain decoupling. 

In Fig. 7(a), we present the staggered AF correlation function between A-spins as J ^ 1. As 
observed, its behavior is well described by that found in a single linear chain, which is asymptotically 
given by 

(-1)' 

C{1) ~ (31) 

apart from logarithmic corrections [s^ . A similar behavior is observed in Fig. 7(b) for the staggered 
AF correlation between Sj-spins up to J = 0.88, a value beyond which the shape of the curve is 
visibly changed. In order to understand this dramatic behavior, we recall that in a two-leg ladder 
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FIG. 6: Magnetic structure factor Fx{q) for the A spins [(a) and (b)], Nc ~ 32, and for the Bi spins 
[(c) and (d) with i = 1 or 2], Nc = 33, for the indicated vahies of J. (e) Magnetic structure factor for 
the composition S; = Bi; + Ba; and J ^ 0.44, 0.45, 0.48, 0.52, 0.56, 0.60, 0.64, 0.68, 0.80, 0.92, 0.96 and 1.00, 
from top to bottom '<3.t q ~ Q. (f) The value of the wave-vector for which the peak at the magnetic structure 
factor exhibited in (e) is observed. Dashed hues are guides to the eye. 



system the asymptotic form of the correlation is given by 



351 



C{1) 



(32) 



where 3.2, see Ref. 



361 ]) defines the correlation length, associated with the gapped spin liquid 
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FIG. 7: Staggered correlation functions between (a) A spins, Nc — 32. and (b) Bi (with i = 1 or 2) 
spins, Nc = 33, for the indicated values of J. In (a) and (b) solid lines indicate the asymptotic be- 
havior for a single chain, (c) Staggered correlation functions between Bt (with i = 1 or 2) spins for 
J = 0.88,0.92,0.96,1.0,1.1,1.2,1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9 and 2.0, circles-dash from top to bottom and 
Nc ~ 33. For comparison we plot the behavior for a single chain (solid line) and for a two-leg ladder with 
32 rungs, (d) Correlation length ^ as a fimction of J: solid line indicates the fitting of the data for Nc = 65 
to Eq. ([55)) with g given by Eq. ([3D]); dashed line indicates the value of ^ (w 3.2) for a two-leg ladder. 



state of this system. Indeed, as displayed in Fig. 7(c) the staggered correlations asymptotically 
approaches the correlation in a two-leg ladder system. In Fig. 7(d) we present the behavior of ^ 
as a function of J for Nc = 33 and Nc = 65. These data were obtained by a proper fitting of Cb^ 
in the interval Iq < I < {Nc/2): starting from J = 2 and taking Iq k, Q (about twice the value of 
^ of a two- leg ladder), we find ^; twice this value of ^ was used as input {Iq = 2^) for the next 
chosen value of J, and so on. Moreover, we have obtained a good fitting to these data by using the 
two-loop analytic form of the 0(3) non-linear sigma model (NLSM) correlation length in (H-l) 



dimension 



3^: 



2w 

^ = ae 3 



9 ) 



(33) 



where a is a constant and g is the NLSM coupling. Further, we assume (see below) that the 
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coupling g is the one suitable to the anisotropic quantum Heisenberg two-leg ladder to the NLSM 



9 = 2k^1 + A (34) 

where J± (Jy) is the exchange coupling between spins at the same rung (leg) and k is a constant 
that depends on the choice of the lattice regularization. 

In order to justify Eq. (j34|) for g, we consider a mapping of the model Hamiltonian, Eq. ([T]), to 
the Hamiltonian of an isolated two-leg ladder by eliminating the spin degrees of freedom associated 
with the A sites. The mapping is performed, in a semiclassical manner, by the following assumption 
on Hab [Eq. 



I 

where 7 is an effective coupling constant. This amounts to reduce the A-B coupling to spins within 
the same unit cell, and cell-cell interactions are taken into account through the effective coupling 
7. We now write: (A; + S;)^ = Af + Sf + 2A; • S;, with Sf = Bfi + B| + 2Bi/ • B2r, since within 
a unit cell (A; + S;)^ (1/2)^ in an AF phase, and dropping constant terms, Hab can be written 
as 

Hab = -1^Bii-B2i. (36) 

I 

Since correlations between spins at A sites does not play a significant role close to the transition, 
we discard the term Ha, and, finally, obtain the following anisotropic two-leg ladder Hamiltonian: 

H^HaL = Hab + Hb, (37) 
where the exchange couplings are given by 

J± = J - 7, (38) 

J|l = J. (39) 
Substituting Eqs. ([MI) and ([39]) into Eq. ([Ml), we find the effective NLSM couphng: 



/ 6(J- Jc2) 

where Jc2 = 7/3. We have fitted the data in Fig. 7(d) to Eq. ([33j) . with g given by Eq. (00]), and 
a, n and Jc2 as fitting parameters. The obtained value of a (=2.7) is such that ^ — > 3.1 as J — > 00, 
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which agrees with the expected value for an isolated isotropic two-leg ladder (~ 3.2), while ac = 4.5 
and Jc2 = 0.91, in agreement with the correlation function behavior shown in Fig. 7(b). 

Finally, in Fig. [8]we display the very interesting behavior of the density of singlets, r], as function 
of J. It is clear that the effect of the A-spins and singlet-singlet interaction is relevant only for 
Jt ^ J ^ Jc2, otherwise the solution, Eq. (21), for low density of singlets can be extended to the 
region of low density of triplets above Jc2 (strongly coupling limit), where correlations between 
B-spins are exponentially small [see Eq. (32)]. In fact, the asymptotic value predicted by Eq. (21), 
i. e., r] = ^/G/tt 0.78, compares well with the numerical one: ~ 0.71. 



V. SUMMARY AND CONCLUSIONS 



Ferrimagnetism 

+ Spiral 
Critical + 

Lieb-Mattis Antiferromagnetic Antiferromagnetic Ladder-Chain 
Ferrimagnetism Transverse Correlations Correlations Decoupling 
I • O • 



Jcl Jf I 



c2 



FIG. 9: Schematic representation of the phase diagram. 



In this work we have derived the rich phase diagram of a three-leg spin Hamiltonian related 
to quasi-one-dimensional ferrimagnets, as function of a frustration parameter J which destabilizes 
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the ferrimagnetic phase. In Fig. [9] we present an ihustration of the obtained phase diagram, 
which displays two critical points, Jd ~ 0.342 and Jc2 ~ 0.91, and a first order transition point 
at Jt ~ 0.445. Through DMRG, exact diagonalization and a hard-core boson model, we have 
characterized the transition at Jd as an insulator-superfluid transition of magnons (built from 
the coherent superposition of singlet and triplet states between B sites at lattice unit cells), with 
a well defined Tonks-Girardeau limit in the high diluted regime. Ferrimagnetism with critical 
staggered correlations in a direction transverse to the spontaneous magnetization is observed for 
Jci < J < Jc2- Further, for Jd < J < Jt the number of singlets in the lattice is quantized, 
while above the first order transition at J = Jt this quantity is a continuous one. Also, in the 
interval Jt < J < Jc2 the magnetic structure factor displays a singlet phase with incommensurate 
(g 7^ and vr) spiral and AF peaks. However, the spiral peak broads and the AF peak is the salient 
feature as J increases within this phase. At J = Jc2 a remarkable gapped two-leg ladder / critical 
single-linear chain decoupling transition occurs, characterized by an essential singularity in the 
correlation length as predicted by the NLSM through a mapping of our model onto an anisotropic 
quantum Heisenberg two-leg ladder. For J ^ Jc2 the ladder approaches the isotropic limit (full 
decoupling), while the linear chain remains critical. 

In summary, our reported results clearly reveal that frustrated quasi-one-dimensional magnets 
are quite remarkable systems to study magnon condensation, including the crossover to coupled 



ladder systems of higher dimensionality 



39j] and related challenging phenomena [4^, as well as 



frustration-driven quantum decoupling transition in ladder systems. 
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